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Abstract 



We present the results of the first lattice study of semileptonic decays of 
baryons containing a 6-quark. Predictions for the decay distributions arc 
given and the Isgur-Wisc functions for heavy baryons are computed, for 
values of the velocity transfer up to about to = 1.2. The computations are 
performed on a 24^ x 48 lattice at /? = 6.2 using the Sheikholeslami-Wohlert 
action in the quenched approximation. 
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1 Introduction 



The discovery of the baryon at LEP and the observation of its semileptonic 
decay 0| makes a study of the weak interactions of heavy baryons on the lattice 
timely. A knowledge of the strong interaction effects in semileptonic decays is 
necessary for the determination of the Vcb element of the CKM matrix from 
the experimentally measured rates and distributions. Up to now \Vcb\ has been 
measured from the inclusive and exclusive decays of heavy mesons. 

In this paper we present the results of the first non-perturbative computation of 
the semileptonic decays of heavy baryons performed using lattice QCD, encour- 
aged by our previous results on the spectroscopy 0. The main purpose of this 
first study is to establish whether such a calculation is feasible, and to identify the 
principal sources of systematic uncertainties and statistical fluctuations. Never- 
theless, in spite of the exploratory character of this investigation, we are able to 
determine the main features exhibited by the six form factors which enter in the 
decay amplitude for the processes A;, — > Ac + W and Ej, —>■ Ec + W, and hence to 
extract a considerable amount of phenomenologically interesting information. 

Using Heavy Quark Symmetry it is possible to relate the six form factors to a 
unique baryonic Isgur-Wise function, ^(cj) Q], where u = v ■ v', and v and v' 
are the four- velocities of the initial and final state baryons. We compute ^{uj), 
and by studying the dependence of the form factors on the mass (or masses) of 
the heavy quark, mg, we are able to use Heavy Quark Effective Theory (HQET) 
to give an estimate of the 0{l/mQ) corrections to the infinite-mass results (for 
a review and references to the original literature see Ref. [^). The dependence 
of the Isgur-Wise function on the masses of the light quarks is more uncertain 
because of our limited data. In particular, we performed the simulations at two 
values of the light-quark mass (a little larger and a little smaller than that of 
the strange quark). Although the results for these masses are very encouraging, 
a significant uncertainty is introduced when the results are extrapolated, as a 
function of the masses of the light quarks, to the chiral limit. We present the 
results for the Isgur-Wise function for massless light quarks (relevant for the decay 
Ab —>■ Ac + W), as well as for the case in which one light quark is massless and 
other is the strange quark, relevant for the process Eb —>■ E^ + Iv. One of the 
main goals of a future simulation will be to determine, in detail, the dependence 
of ^(cu) on the masses of the light quarks, and hence to reduce the uncertainty 
due to the extrapolation of the results to physical masses. 

We obtain for the decay rate integrated over the range uj G [1, 1.2] 

1-1.2 

du—iAb^Ac + W) = lA tl\Vcb\' 10''' 

du—{Eb^Ec + lv) = l.Q tt\Vcb\'w''s-\ (1) 
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This range of uj corresponds to that for which we have the most rehable results. 
We also obtain the slopes of the differential-decay-rate form factors which are to 
be compared to the slopes obtained by performing fits to experimental results for 
dV/duj versus a;, for a; near 1, when such results become available. We find: 

p| = l.l±1.0 (2) 

for Af, kciv decays and 

P0 = 1.4±O.8 (3) 

for Sft — > 'Ef.iv decays. A more detailed discussion is presented in Section ^]2[ 

The plan of the paper is the following: in Section 2 we recapitulate the theoret- 
ical framework for the definition of the baryonic Isgur-Wise function. Section 3 
gives details of the simulation and Section 4 presents the results of the numerical 
analysis of both two- and three-point functions. We study the dependence of the 
Isgur-Wise function on the velocity transfer and on the masses of the heavy and 
light quarks in Section 5. In Section 6, we consider the phenomenological im- 
plications of our results for the baryonic Isgur-Wise function and give estimates 
of differential and partially-integrated decay rates. The main body of the paper 
is accompanied by three appendices where we lay out some lengthy parts of the 
calculation and of the analysis. 

For the reader who is not interested in the details of the computation we have 
attempted to write Sections 2 (theoretical framework) and 6 (results and impli- 
cations for phenomenology) in a self-contained way. 

2 Theoretical background 

In this section we outline the theoretical framework needed to study the semilep- 
tonic decays of A^ and Sf, baryons. The two baryons differ in their light flavour 
content, but are identical in all other quantum numbers and, in particular, they 
can be described in the same manner within the framework of HQET. For sim- 
plicity we will only consider the A;, decay. 

The non-perturbative strong interaction effects in the exclusive semileptonic de- 
cay of the Afc are contained in the matrix elements of the V — A weak current, 
Jfi = C7^(l — 75)^, which can be written in terms of six invariant form factors, 
fi, Qi with i = 1,2, 3, as follows: 

{Aif\p')\J,\At\p)) = u[^\p') Yl,Ui - l^gi) + ^cy,,q^{h - 75^72) 

+^q,{h-l,g,)]ut\p). (4) 

In Eq. (^, the momenta and polarisations of the initial and final baryons have 
been explicitly indicated. The form factors are functions of q^, where q is the 
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4- momentum transfer {q = p' — p)- The decomposition above is convenient since, 
if the lepton masses are neglected, only the dominant form factors /i and gi 
contribute to the rate. Since both the quarks destroyed and created by are 
heavy, HQET provides a useful guide to the study of the form factors. In addition, 
the Afc baryon has a particularly simple structure in that it is composed of a heavy 
quark and light degrees of freedom with zero total angular momentum, so that 
Heavy Quark Symmetry has considerable predictive power. Expression can 
be rewritten in terms of the velocity variables = p^/M\^ and v' = p'/M, 



Ac 



u = V ■ v': 



using a different set of form factors, which are functions of the velocity transfer 



4'\v') [7^(Fi - 756-1) + v^iF2 - 7562) + v'^iFs - 75G3)] ut\v). (5) 

In principle, we do not need to use HQET at all but can evaluate the matrix ele- 
ments in Eq. (^) directly. However, as the inverse lattice spacing in our calculation 
is about 2.9 GeV, we cannot use the physical 6-quark mass in our computation. 
Instead, we evaluate the matrix element in Eq. (^) for a variety of "6" and "c" 
quark masses in the region of the physical charm quark mass, and extrapolate the 
results to the physical values. HQET provides us with the theoretical framework 
to perform this extrapolation. The hadronic form factors Fi, Gi can be expanded 
in inverse powers of the heavy-quark masses; the non-perturbative QCD effects 
in the coefficients of this expansion are universal, mass-independent functions 
of the velocity transfer. This is achieved by constructing an operator product 
expansion of the weak currents. The analysis at leading order establishes 
the important result that all the baryonic form factors are described by a single 
universal function called the (baryonic) Isgur-Wise function, ^{to, fj.), 

(A(^)(t;')|c7.(l - l6)b\At\v)) = M'\v'hua ' l^>t\v), (6) 

which is normalized to the identity at = 1 and where fi is the scale at which ^ 
is renormalized. The relation between the six form factors and the renormalized 
Isgur-Wise function is given at leading order by the expressions p, ^ 

Fiioo) = Q(.;)f-(c.), G,,{uj) = (c.)f-(u;), (7) 

where the scale dependence of the Isgur-Wise function ^{uj,fi) is reabsorbed into 
the definition of the short distance coefficients Ci^\uj). The coefficient functions 
C*j-^^ are known up to order 0^(2; log 2)", where z = mc/mi, is the ratio of the 



^We draw the reader's attention to the fact that the Isgur-Wise function which describes 
baryonic weak matrix elements is different from that entering in mesonic transitions. For 
simphcity we shall use the name Isgur-Wise function to refer to the baryonic one, whenever 
this does not lead to ambiguities. 
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heavy-quark masses and n = 0, 1, 2. The perturbative expansion of Ci and Cf is 
of the form 

Ci"^ = l + 6fiuj)as + ... (8) 

whereas the expansion of Cg^j starts at 0{as)- The perturbative corrections are 
rather small since the coefficients of as are typically of order 1, for our range of 
masses and values of tu, so that C2 3 are much smaller than (Cfj ~ 0.2 (7^). 

We also attempt below to obtain some information on the l/mq corrections to 
the form factors. In order to include these corrections it is convenient to define a 
new function (see Ref. M) 



A 



A 



2m', 



Q 



UJ 



[UJ 



(9) 



which is also scale-independent and normalized at zero recoil, since the func- 
tion x{^)y arising from the higher- dimension operators in the HQET Lagrangian, 
vanishes at u = 1. 



A in Eq. 

baryon, 



is the binding energy of the heavy quark in the corresponding A- 



A = Ma. -mb = M^, - nic 



(10) 



up to l/rriQ corrections. The quark masses in Eq. (|I0|) are generally taken to be 
pole masses, which contain renormalon ambiguities of 0{Aqcd)- For the form 
factors, the ambiguities due to A are cancelled by those arising from the higher- 
order terms in the perturbative series for the coefficients {00). In practice, 
we only know the coefficients Ci^\uj) up to one-loop order, and we implicitly 
assume that A is obtained from some physical quantity with a similar precision. 
As we note in Section ^ some of the freedom of what to assign to the coefficients 
and what to the power corrections cancels in the prediction for the physical form 
factors obtained from those calculated on the lattice. 

^QQ' is not a universal function and its dependence on the flavour of the heavy 
quarks must be studied in detail (see Section . The relation between the new 
function C,qq' (uj) and the form factors is given by 



F,{uj) = Ni{uj)^QQ,{uj) + 0[l/m'Q,,, 



where the coefficients N-^^ contain both radiative and l/mq corrections. The 
exact expressions are p[ 



iVi(^) 



2 / A A 



a; + 1^2mQ 2mQ'\ 



5 



N2{UJ) 



C2(CU) (^1 + 

Csioo) (l + 



2uj A 
u + 1 2mQ 

2uj A 
oj + 1 2mr 



Ci(cu)+C3(c^) 



A 



uj + 1 2mQ 
2 A 



LU + 1 2m, 



Q 



'12) 



(c^) 1 + 



2 A A ^ 

uj+12mQi mq^ 



A 



uo + 1 2mr 



A 



+ 



A 



to +1 2mQ niQf J 
2 A 
u + 1 2mQ 

where the "velocity transfer of the free quark" uj is given by 



Cf (^) - 



, / A A 

UJ = UJ + [uj — 1) \ h 



mq TUq, 



(13) 



At zero recoil, Luke's theorem [0 protects the quantities J2i -^i(l) and G'i(l) from 
0{l/mQ) corrections: 



rrir 



G,{l)=C!{l) + 0{l/ml) 



iVsum(l) 

Gijl) 
iVf(l) " 



= l + 0(l/m^) 

l + 0(l/m^). (14) 



For degenerate quark masses, tuq = rriQ/, the coefficient 



N,U^) = Y,N,{uj) = Ci{uj) + C2{uj)[l 

i \ 
+ C^{oj){l 



2{uj-l) A 

UJ + \ 2mQ ^ 
2iuj-l) A 



(15) 



a; + 1 2mQi J 

equals one at zero recoil, as required by vector current conservation. 

The coefficients (|T2|) and (|T3|) depend on the heavy-quark masses and on the value 
of the baryonic binding energy A. Different choices for the definition and values 
of the quark masses to be used in these expressions lead to differences which are 
of 0{al) and 0{l/mQ), and hence are formally of the size of the terms we are 
neglecting. Moreover, we find that different choices of the quark mass lead to 
negligible differences in the form factors for Fi, and Gi. This is not the 
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ivip a 




IVlv CL 




mQ[(je\ ) 


0.121 


0.874 1 


4 
3 


0.896 1 


5 
4 


2.38 


0.125 


0.773 1 


3 
3 


0.799 :^ 


4 
3 


2.10 


0.129 


0.665 1 


3 
3 


0.696 :^ 


4 
4 


1.80 


0.133 


0.547 1 


3 
3 


0.588 t 


4 
5 


1.48 



Table 1: Values of the quark masses in physical units, for = 2.9 GeV and 
A = 200 MeV, as obtained from the pseudoscalar and vector meson masses in 
lattice units, (kq is the heavy quark's hopping parameter.) 

case for the form factors -^2,^3 and 02,0^, however, for which the coefficient 
functions are zero at tree level and for which the 0{l/mQ) terms represent a 
major contribution. In these cases we take mq = M\ — A, and A is obtained 
from a fit of the theoretical prediction for the form factors to the lattice results 
(see Section ^.3|) . In all other cases, where the results are insensitive to the choice 
of the quark mass, we take 

mg = — (3My + Mp) - 200 MeV (16) 

where My and Mp are the masses of the vector and pseudoscalar heavy-light 
mesons, as measured very precisely from a previous simulation performed on the 
same set of configurations. 200 MeV is an estimate of the spin-averaged mesonic 
binding energy 0. For the lattice spacing, as discussed in greater detail in Section 
^, we take 

= 2.9 ± 0.2 GeV. (17) 

The values of the quark masses obtained in this way are presented in Table |I|, and 
the coefficients N^^^ and C*, which are used in the rest of this paper are reported 
in Tables ^ and |^ for Aqcd = 250 MeV and ri/ = (as may be more appropriate 
for a quenched calculation). 

Below, we will proceed as follows. Since Nf and A^^sum are close to unity and the 
form factors Gi{uj) and J^i^^iiuj) are insensitive to the value of A, we can use 
the measured values of these form factors to determine the Isgur-Wise function 
reliably. For the individual vector form factors Fi, F2 and F3 (and similarly for G2 
and G3) the situation is different since the 0{l/mQ) corrections are significant. 
From these form factors, using the Isgur-Wise function already obtained, we 
determine the A parameter. Finally we use the combined results to determine 
the form factors for the physical b ^ c decays. 
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0.121 
0.121 
0.121 

0.121 



0.121 
0.125 
0.129 

0.133 



1.0 
1.0 
1.0 
1.0 



1.0 
1.0 
1.0 
1.0 



0.961 
0.965 
0.971 

0.98 



0.019 
0.024 
0.03 

0.04 



0.961 
0.965 
0.971 

0.98 



1.0 
1.01 
1.02 

1.04 



0.125 
0.125 
0.125 
0.125 



0.121 
0.125 
0.129 
0.133 



1.0 
1.0 
1.0 
1.0 



1.0 
1.0 
1.0 
1.0 



0.956 
0.959 
0.965 
0.973 



0.016 
0.02 
0.027 
0.037 



0.956 
0.959 
0.965 
0.973 



0.992 
1.0 
1.01 
1.03 



0.129 
0.129 
0.129 
0.129 



0.121 
0.125 
0.129 
0.133 



1.0 
1.0 
1.0 
1.0 



1.0 
1.0 
1.0 
1.0 



0.949 
0.952 
0.957 
0.964 



0.011 
0.016 
0.022 
0.031 



0.949 
0.952 
0.957 
0.964 



0.982 
0.989 
1.0 
1.02 



0.133 
0.133 
0.133 
0.133 



0.121 
0.125 
0.129 
0.133 



1.0 
1.0 
1.0 
1.0 



1.0 
1.0 
1.0 
1.0 



0.938 
0.941 

0.945 
0.952 



0.0044 
0.0086 
0.015 
0.024 



0.938 
0.941 

0.945 
0.952 



0.967 
0.973 

0.984 
1.0 



0.121 
0.121 
0.121 

0.121 



0.121 
0.125 
0.129 

0.133 



1.1 
1.1 
1.1 
1.1 



1.15 
1.15 
1.16 
1.17 



0.932 
0.937 
0.944 

0.955 



0.018 
0.022 
0.028 
0.037 



0.932 
0.937 
0.944 

0.955 



0.968 
0.977 
0.99 
1.01 



0.125 
0.125 
0.125 

0.125 



0.121 
0.125 
0.129 

0.133 



1.1 
1.1 
1.1 
1.1 



1.16 
1.16 
1.17 



0.927 
0.931 
0.938 

0.918 



0.014 
0.018 
0.024 



0.927 
0.931 
0.938 

0.918 



0.96 
0.969 
0.981 

1.0 



0.129 
0.129 
0.129 
0.129 



0.121 
0.125 
0.129 
0.133 



1.1 
1.1 
1.1 
1.1 



1.16 
1.16 
1.17 
1.17 



0.92 
0.924 
0.93 
0.94 



0.01 
0.014 
0.02 
0.028 



0.92 
0.924 
0.93 
0.94 



0.95 
0.958 
0.97 
0.99 



0.133 

0.133 
0.133 
0.133 



0.121 

0.125 
0.129 
0.133 



1.1 
1.1 
1.1 
1.1 



1.17 
1.17 
1.17 
1.18 



0.91 

0.913 
0.919 
0.929 



0.0039 

0.0077 
0.013 
0.021 



0.91 

0.913 
0.919 
0.929 



0.935 

0.943 
0.954 
0.973 



Table 2: Radiative and 0{l/mQ) correction factors for uj = 1.0 and 1.1. (kq and 
kqi are the initial and final heavy quark hopping parameters.) 
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UJ 



UJ 



0.121 
0.121 
0.121 

0.121 



0.121 
0.125 
0.129 

0.133 



1.2 
1.2 
1.2 
1.2 



1.3 
1.31 
1.32 
1.33 



0.904 
0.909 
0.917 

0.929 



0.016 
0.02 
0.025 

0.034 



0.904 
0.909 
0.917 

0.929 



0.937 
0.946 
0.96 

0.981 



0.125 
0.125 
0.125 
0.125 



0.121 
0.125 
0.129 
0.133 



1.2 
1.2 
1.2 
1.2 



1.31 
1.31 
1.32 
1.34 



0.898 
0.903 
0.911 
0.923 



0.013 
0.017 
0.022 
0.03 



0.898 
0.903 
0.911 
0.923 



0.929 
0.938 
0.951 
0.972 



0.129 
0.129 
0.129 
0.129 



0.121 
0.125 
0.129 
0.133 



1.2 
1.2 
1.2 
1.2 



1.32 
1.32 
1.33 
1.35 



0.891 
0.896 
0.903 
0.915 



0.0091 
0.013 
0.018 
0.026 



0.891 
0.896 
0.903 
0.915 



0.919 

0.928 
0.941 
0.961 



0.133 
0.133 
0.133 
0.133 



0.121 
0.125 
0.129 
0.133 



1.2 
1.2 
1.2 
1.2 



1.33 
1.34 
1.35 
1.36 



0.881 
0.886 

0.893 
0.904 



0.0034 
0.0069 

0.012 
0.019 



0.881 
0.886 

0.893 
0.904 



0.905 
0.913 

0.925 
0.945 



0.121 
0.121 
0.121 

0.121 



0.121 
0.125 
0.129 

0.133 



1.3 
1.3 
1.3 
1.3 



1.45 
1.46 
1.48 
1.5 



0.876 
0.882 
0.89 

0.903 



0.015 
0.018 
0.023 

0.031 



0.876 
0.882 
0.89 
0.903 



0.907 
0.916 
0.93 

0.951 



0.125 
0.125 
0.125 

0.125 



0.121 
0.125 
0.129 

0.133 



1.3 
1.3 
1.3 

1.3 



1.46 
1.47 
1.49 
1.51 



0.871 
0.876 
0.884 

0.897 



0.012 
0.015 
0.02 

0.028 



0.871 
0.876 
0.884 

0.897 



0.899 
0.909 
0.922 

0.913 



0.129 
0.129 
0.129 
0.129 



0.121 
0.125 
0.129 
0.133 



1.3 
1.3 
1.3 
1.3 



1.48 
1.49 
1.5 
1.52 



0.864 
0.869 
0.877 
0.889 



0.0083 
0.012 
0.016 
0.024 



0.864 
0.869 
0.877 
0.889 



0.89 
0.898 
0.911 
0.932 



0.133 

0.133 
0.133 
0.133 



0.121 

0.125 
0.129 
0.133 



1.3 
1.3 
1.3 
1.3 



1.5 
1.51 
1.52 
1.54 



0.854 

0.859 
0.866 
0.878 



0.003 

0.0062 
0.011 
0.018 



0.854 

0.859 
0.866 
0.878 



0.875 

0.884 
0.896 
0.916 



Table 3: Radiative and 0{l/mQ) correction factors for uo = 1.2 and 1.3. (kq and 
kqi are the initial and final heavy quark hopping parameters.) 
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3 Details of the simulation 



Our calculation is performed using 60 SU (3) gauge field configurations generated 
on a 24^ x 48 lattice at /3 = 6.2, using the hybrid over-relaxed algorithm described 
in Ref. 0. Since we are studying the decays of quarks whose masses are large in 
lattice units, we must control discretisation errors. In order to reduce these errors, 
we use an O (a)-improved fermion action originally proposed by Sheikholeslami 
and Wohlert|TO|, given by 

= S^-t^Y. q{x)F,,{x)a,,q{x), (18) 
where Sp' is the Wilson action: 

= ^Ui^)'li^) - «^Il[^(a;)(l -7M)f^M(a;)g(a; + /i) 

X \ fj. 

+ g(x + /i)(l + 7^)t/t(x)g(x)]j . (19) 

The use of the SW action reduces discretisation errors from 0{ma) to (9(asma) 
|I0| , |1^] provided one also uses "improved" operators, for example, those obtained 
by "rotating" the field of the heavy quark, Q: 

Q(x)^il-^rD)Q{x) . (20) 

Thus, to obtain an O (a)-improved evaluation of the matrix element in Eq. (^, 
we use a "rotated" improved current 

= Q'{x)t^Q{x) , (21) 

where 

f = (1 + It D) (1 - 7^^)75 (1 - \l- D). (22) 

The gauge field configurations and the light quark propagators were generated on 
the 64-node i860 Meiko Computing Surface at the University of Edinburgh. The 
heavy quark propagators were computed using the Cray T3D, also at Edinburgh. 

In order to enhance the signal for the baryon correlation functions, the light and 
heavy quark propagators have been computed using the Jacobi smearing method 
[0, at both the sink and the source (SS). Since smearing is not a Lorentz- 
invariant operation, it alters some of the transformation properties of computed 
quantities. In a previous publication 0, we have shown that such an effect is 
evident in two-point baryonic correlators at non-zero momentum. In Appendix 
0, we include a study of the smearing effects for SS three-point functions. 
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Statistical errors are obtained from a bootstrap procedure [|13|. This involves the 



creation of 1000 bootstrap samples from the original set of 60 configurations by 
randomly selecting 60 configurations per sample (with replacement). Statistical 
errors are then obtained from the central 68% of the corresponding bootstrap 
distributions as detailed in Ref. [|T3 



3.1 Correlators and details of the analysis 

In order to study semileptonic decays of the type Af, Ac / ui on the lattice, we 
consider the following three-point correlators: 

iC{t.,ty)f,-'^' = 5:^e-^"-V^«^-^^(O«'(x)(J^(i/))«-«'(5«(0)) (23) 

X y 

where the spinorial indices are implicit. The operator 

0'3ix)=eaUlfCl5ll)Q'' (24) 

is the interpolating operator for the A baryon, and the current which mediates 
the decay of a heavy quark Q into a second heavy quark Q' is given in Eq. (^Tj). 
h and I2 represent light quark fields. The three-point function {C{tx,ty))^^'^ 
can be written in terms of quark propagators, as 

(C(t.,t,))«-«' = -{Y^{E^'\0,y;t,,p)Q,T,S$\y,0))e-^'^-y^, (25) 

y 

where Sq^y, 0) is the propagator of the Q quark from y to the origin in the 
presence of a background field configuration, 

E^'^(0,i/;t.,p)Q, = e,,,e,,,,,,^e-^P"-^^Tr[^f ^(x, 0)C75^f (x, 0)75^] ^^'^(x, y) 

X 

(26) 

where T represents the transpose in spinor-space, (..) denotes the average over 
gluon configurations, and a, b, c, d, a', b', d are colour indices. 5*/^ and Si^ are 
the propagators of the two light quarks. The extended propagator (120) can be 
evaluated using the standard source method reviewed in Ref. |jl5|. 

In the limit of large tx and — ty (in the forward part of the lattice), where the 
ground state contributions to the correlation function should dominate, we can 
rewrite the correlator as follows i 

7 7' 

iC{tx, ty)fr''' = {i^' + ^0^?' ^(I^ + M)) , (27) 



^This and the following expressions are only correct in the case in which local operators are 
used. The actual case of smeared-smeared correlators, which is discussed in Appendix ^ is 
more complicated but conceptually similar. 
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where 

E = + \^ E' 

M = Ma(q) M' 

Z = Z{A{Q),\p\) Z' 
The weak matrix element can be written in terms of the six lattice form factors, 

(Ag(p 01(^.(0, 0))«-«'|A«(p)) = ^S'Hp')-^?''^(p',p)«S3^Hp) (29) 

with ^f'^(p',p) = {fI{uj)i^^F^{u)v^^FI{uj)v'^ 

- (GfM7, + G^Mt;, + G^Ht;;) 75 (30) 

is] frl 

where ^q(q/) are the spinors of the heavy baryons and s and r are helicity indices. 
The lattice form factors and G\ are related to the physical form factors Fj and 
Gi through the renormalisation constants of the lattice vector and axial current, 
respectively. This will be discussed in some detail below. 

4 Details of the numerical analysis 

In this section we describe our procedure for the extraction of the form factors 
from the lattice correlation functions. 

4.1 Analysis of two-point correlation functions at non- 
zero momentum 

The wave-function factors Z, Z' and the energies which appear in Eq. ( pTj) can 
be obtained from the analysis of the appropriate two-point correlation functions: 

CiiM = Ee-^^^-^^(O«(f,t)O'^(0,0)). (31) 

This correlator was evaluated for four values of the heavy-quark hopping pa- 
rameter corresponding to masses around that of the charm quark, and for three 
combinations of the light-quark masses, as shown in Table (^ . We have computed 
the correlator for momenta up to \p\ = 2pniin, where Pmm = ^n/La = n /Via is the 
minimum non-zero momentum allowed on our lattice. For the analysis, equivalent 
momenta have been averaged to reduce the statistical fluctuations. 

For the actual case of smeared-smeared operators (see Appendix 0), for large time 
separations and imposing antiperiodic boundary conditions in time, the correlator 



= Ma(q') /28 
= Z{A{Q'),\p-q\). 
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0.14144/0.14144 


0.121 


0.125 0.129 0.133 


0.14144/0.14226 




- 0.129 - 


0.14226/0.14226 




- 0.129 - 



Table 4: Quark hopping parameter combinations used in the calculation of baryon 
two-point functions, (kh and k,i2 are the two light-quark hopping parameters and 
kq, that of the heavy quark.) 



C9 b 



ecomes 



zm){ 



C^{t,P} 



e ^* 



E+M-a^{E-M) ^ . E+M+a^{E-M) _ r; 



-E(T-t) 



\ E+M-a^{E-^M) ^ E+M+aHE~M)_ 2a- -1 \ 



where p^ is defined in (p8|), p^ = {E, —p) is the 4- momentum of the antibaryon 
propagating in the backward part of the lattice, and Zg, a are the amplitudes of 
the smeared operator (defined in Eqs. (|8ll) and (^) ). 

At zero momentum, the analysis is particularly simple since the only non-zero 
components are independent of the unphysical amplitude a. Details of this case 
are given in Ref. f^. At finite momentum, we have used both the diagonal and 
the off-diagonal components of the spinorial matrix to extract the energy from 
the exponential fall-off. The values of the amplitudes Zg and a were obtained by 
fitting separately the contributions proportional to the identity, to 70 and to 7^, 
and by taking suitable linear combinations of the overall factors, as explained in 
the appendix of Ref. |Q. The results of these fits, for those mass combinations 
which are relevant to the present study, and for momenta up to \p\ = \^n/12a 
are reported in Tables |^ and ^. 

The case corresponding to kh = = 0.14144 is further studied in detail to check 
the precision with which the dispersion relations are satisfied as the momentum 
of the baryon is increased. It is nowadays customary [0, to replace the 
continuum dispersion relation (cdr) 

a'^E^ = a'^m'^ +p'^a'^ (33) 

with the so-called lattice dispersion relation (Idr) 

a^E^ = a^m^ + sin^ (pa), (34) 

which is suggested by the form of the free fermionic propagator on a discrete 
lattice. For the heavy-quark masses in Table ^ the two dispersion relations yield 
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p= (0,0,0) 


P - li^' 


n f\\ 
u, UJ 


P - ll^' 


2tt n\ 


14144/0 14144 


0.121 


1.138 


+ 7 

T 
— 1 


1.167 


+ 9 

-9 


1.192 


+ 13 
-12 




0.125 


1.040 


+ 6 
-6 


1.072 


+ 9 
-9 


1.099 


+ 13 

-12 




0.129 


0.938 


+ 6 
-6 


0.973 


+ 8 


1.002 


+ 13 




0.133 


0.829 


+ 6 
-6 


0.868 


+ 8 
-7 


0.901 


+ 14 
-13 


0.14144/0.14226 


0.129 


0.910 


+ 8 
-7 


0.943 


+ 9 
-9 


0.972 


+ 12 
-14 


0.14226/0.14226 


0.129 


0.876 


+ 9 
-8 


0.914 


+ 10 
-10 


0.948 


+ 17 
-17 


chiral / chiral 


0.129 


0.807 


+ 15 
-10 










chiral/strange 


0.129 


0.853 


+ 14 
- 6 











Table 5: Energies of the A-baryon in lattice units for all the momenta and quark 
hopping parameters relevant to the present study. 





p= (0,0,0) 


P= (Pmin,0,0) 


P= (Pmin,Pmin,0) 




X 10^ 


Z'^ X 10* a 


Z'i X 10* a 


{h-n. h-i:2) = (0.1 11 IL 0.1 11 11) 


0.121 
0.125 
0.129 

0.133 


4.44 ttl 
4.41 tt^ 
4.35 ttl 

1.17 -^^ 


2.85 ttl 0.66 

2.86 0.70 
2.84 til 0.77 tl 

2.75 t^ 0.83 tl 


1.96 ttl 0.56 1^ 
1.95 tt' 0.60 
1.94 0.66 t^l 

1.87 0.73 +1? 


{nii,fii2) = (0.14144,0.14226) 


0.129 


4.02 ttl 


2.56 til 0-76 tl 


1.68 0.68 1}? 


(«:n,«:/2) = (0.14226,0.14226) 


0.129 


3.77 1:% 


2.42 til 0.71 tl 


1.60 ttl 0.65 1}^ 



Table 6: Amplitudes Zg and a obtained from the analysis of the finite momentum 
two-point functions. 
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essentially indistinguishable results for the energy at |p| = Pmin- In addition, 
the theoretical predictions coincide with the measured values, confirming that 
the systematic effects at this low value of the momentum are negligible. For 
momentum \p\ = y^pmin, we note that the predicted value is about 1% larger 
than the measured one, although always compatible within one sigma. Given this 
result, it seems that the correction obtained with the lattice dispersion relation 
goes in the right direction. However, much more precise data are needed to draw 
a firm conclusion on this issue. 

Finally, for the conversion of our values for masses and energies into physical units 
we need an estimate of the inverse lattice spacing in GeV. Following Ref. [Q, we 
use the value given in Eq. (|T^ . The error in Eq. ( ]T7| ) is large enough to encompass 
all our estimates for from quantities such as rrip, /tt, itln, the string tension 
\/K and the hadronic scale Rq discussed in 



4.2 Three-point functions and lattice form factors 

In this subsection we explain our procedure for extracting the form factors from 
the computed three-point (and two-point) correlation functions. We have com- 
puted the three-point functions for the mass combinations tabulated in Table 0. 
In order to study the dependence of the form factors on the masses of the heavy 
quarks, we have computed the correlation functions for all combinations of kq 
and K,Q' taken from 0.121, 0.125, 0.129, 0.133, but with the light-quark masses 
fixed by kh = K12 = 0.14144 (which is close to that of the strange quark). On the 
other hand, the dependence on the hght-quark masses was studied by keeping 
fixed the mass of the initial and final heavy quarks kq = kq/ = 0.129 which is 
very close to that of the charm quark, and considering the three light hopping- 
parameter combinations kh = K12 = 0.14144, kh = 0.14226,^^2 = 0.14144, and 
ft^zi = 1^12 = 0.14226. In light of the encouraging results obtained below, we envis- 
age the possibility of repeating the calculation on a larger sample of heavy- and 
light-quark masses. 

We measure the six form factors for the different quark masses from the three- 
point functions (p3|), whose expression, for large values of tx and tx — ty in the 
forward part of the lattice and for local interpolating operators was given in 
Eq. ([27|). The more complicated case of SS correlators, which is explained in 
detail in Appendix 0, can be written schematically as follows: 



ZZ' 
16EE' 

+ W(a, AM,M',p,p')(^XM +<i^3^H 
- A^{a,P,M,M',p,p')G^{u;) 

+ B{a,(3,M,M\p,p'){v,G^{uj) + v'^Gf;{uj))] (35) 
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degenerate transitions 


0.14144/ 
0.14144 


0.121^0.121 0.125^ 0.125 0.129 -^ 0.129 0.133 ^0.133 


0.14144/ 
0.14226 


0.129-^0.129 


0.14226/ 
0.14226 


0.129^0.129 


non-degenerate transitions 


0.14144/ 
0.14144 


0.121 -> 0.125 0.121 -^ 0.129 0.121 ^0.133 
0.125^0.121 - 0.125^ 0.129 0.135 -^0.133 
0.129^0.121 0.129^ 0.125 - 0.133 -^0.133 
0.133-^0.121 0.133^ 0.125 0.133 -^0.129 



Table 7: Quark hopping parameter combinations used in the calculation of baryon 
three-point functions, (nn and k,i2 are the two light-quark hopping parameters 
while kq and kq/ are those of the initial and final heavy quarks.) 

where we have set tx = T/2 (which is the value we chose in our computations), 
and have renamed the remaining time variable ty t. The matrices V, W, A, B 
are different linear combinations of the Dirac matrices, with coefficients which 
depend on the arguments shown. The overall normalisation factor, as well as the 
four matrices V, W, A, B can be fully reconstructed from the two-point functions 
at the corresponding values of the masses and momenta, as illustrated in Section 
p. 

We will now consider separately the two cases in which we keep either the axial 
or vector currents from the V — A current which mediates the semileptonic weak 
decays of heavy baryons. The following discussion is based on equations (^3]) and 
(IID in Appendix 0. 

4.2.1 Analysis of the axial form factors 

The freedom to choose the Lorentz index /i and the spinor coniponents appropri- 
ately allows us to extract the required form factors efficiently □. It is convenient 
to think of the four- by-four spinorial matrix as subdivided into four two- by-two 
matrices, as explained in Appendix We find that 



■^Alternatively we could use suitable projection operators in spinor space. 
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1. current indices fi = i with i = 1,2,3: the large componentgj of Ai are 
located in the top (bottom) diagonal submatrix in the forward (backward) 
part of the lattice. The large components of B are located in the top-right 
(bottom-left) corner submatrix. Thus the contributions of the form factor 
Gi can be separated from that of G2 3- 

2. current index /i = 0: the large components of both ^0 and B are located 
in the top-right (bottom-left) corner matrices in the forward (backward) 
part of the lattice. So the three form factors contribute to the same spino- 
rial components, making the extraction of the form factors very uncertain. 
Therefore, the equations obtained from the current with index /i = will 
not be considered in the analysis below. 



We conclude that from an analysis of the correlators with index /x = 1, 2, 3, it is 
possible to obtain a clean straightforward determination of the form factors. We 
consider the asymptotic form of Eq. (^), fitting for either or 3 separately, 
depending on the particular spinorial component under study. 



Once the time-dependent factor in ( P3[ ) has been divided out, we observe long and 
stable plateaux, centered around t = 12. As an example we exhibit in Figures |l]A 
and |I|B, the plateaux for the form factors and G2, for one set of masses and 
choice of momenta. We have fitted for 3, 5 and 7 timeslices. The central 
values of the fits are insensitive to the choice of fitting interval, and in order to 
avoid possible contaminations due to exited states, we have decided to restrict 
the fitting range to three time slices centered around t = 12. The Xdof these 
fits were always very reasonable, ranging from Xdof ~ 0-5 to Xdof ~ 2. 

As was discussed in Section ^, the form factors are related to the physical renor- 
malized Isgur-Wise function through a multiplicative renormalisation which takes 
into account the short-distance QCD corrections. Furthermore an additional 
renormalisation constant (Zj^) must be introduced in order to relate the lat- 
tice improved axial current (^) to the continuum one. Perturbative and non- 



perturbative calculations of Za are available in the literature [0, but it can 
also be estimated non-perturbatively in this computation. As can be seen from 
Eqs. (|14D) Gi reduces to the Isgur-Wise function, at zero recoil, and it is free 
from 0{l/mQ) corrections because of Luke's theorem. If we make the reasonable 
assumption that 0{l/mQ) corrections are small, a measure of G^{l)/Nf{l) will 



^Here and in the following, we refer to the numerical coefficients of the form factors in 
Eq. ( ^5[ ) as large if they are proportional to the energy or mass of the baryon, and as small if 
they are proportional to the spatial momentum p. 
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-0.1 



-0.2 




D 



Figure 1: Examples of the plateaux used in the fit to (from left to right, top to bot- 
tom) Gi{uj), G2{uj), J2iFi{uj) and ^2(0;). nn = = 0.14144, kq = kq' = 0.121, 
the initial particle is at rest and the final particle has momentum {p', 0, 0), where 
p' = 2n/La. 
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give us a non-perturbative estimate of Zj^. For degenerate transitions, 
iVf(l) 



0.98 +8 at = = 0.14144, and Kh = 0.133 

- 0.97 at Kii = Ki2 = 0.14144, and = 0.129 



Gf(l) • ' 

G[(l) 
Nfil 



(36) 



0.98 +8 at Kn = k/2 = 0.14144, and Kh = 0.125 



G[(l) • 
Gf(l) 



0.98 +^[! at K/i = = 0.14144, and kh = 0.121 



to be compared with the non-perturbative estimate of Za, obtained for hght 
quarks pO[ 

^non-pert ^ ^ +1_ ^3-^^ 

The coefficients Nf used in Eqs. (^) were computed from expressions (p!2D. 
Within the statistical precision of our calculation, we have no evidence of dis- 
cretisation errors nor of 0{l/mQ) corrections (assuming that these two effects do 
not partially cancel, which is very unlikely give the range of masses we consider). 
Indeed, we observe that the four values in ( P^D are exceptionally stable with the 
quark mass. 

In Tables ^, ^ and |TD|, we present our results, for the quantity 

M^) - ^ X ^ (38) 

for all the quark masses and for initial and final momenta up to \p'\ = It 
follows from the above discussion that this quantity is independent of the lattice 
renormalisation constant Za- 
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Ka = 0.14144, K,/2 = 0.14144 



[Pmii 



0.121 0.121 



0.121 0.125 



a; 



0.121 ^ 0.129 



0.121 ^ 0.133 



(0,0,0), 
(1,0,0) 

(0,0,0), 

(1,1,0) 

(1.0.0). 

(0,1,0) 

(1,0,0), 

(1,0,0) 

(1,0,0), 
(-1,0,0) 

(1,0,0), 
(0,0,0) 

(1,0,0), 
(0,1,1) 



1.026 0.97 



1.050 0.95 



1.052 0.94 



+ 3 
-3 



+ 6 



-11 
-12 



1.030 0.96 ^3 



1.060 0.93 



1.057 0.92 



1.000 1.02+1? 



1.10 0.81 +o 



1.077 0.95 



+ 2 



+ 7 



-10 
-10 



1.000 0.95+12 



1.11 0.76^8 
1.026 0.95+^ 



1.09 0.92+}^ 



1.037 0.94+3 



1.073 0.91 1^ 



1.064 0.91 



1.001 0.94+}^ 
1.13 0.74 +? 



1.026 0.95+^ 



1.10 0.88+}^ 



1.048 0.91 ^4 

1.09 0.90 1^ 

1.07 0.89 1^° 

1.003 0.93+12 

1.15 0.70 

1.026 0.95+5 



1.12 0.80 1}^ 



Table 8: Estimates of the function ^qq> as obtained from the axial form factor 
Gi. All the transitions, corresponding to initial heavy k = 0.121 and 0.125 and 
initial and final momenta up to \p\, \p'\ = \/2pmin o,re shown. Statistical errors 
in uj are in the last digit or beyond. 
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Kn = 0.14144, Kz2 = 0.14144 



bmir 



0.125 ^ 0.125 



0.125 0.121 



0.125 0.129 



0.125 ^ 0.133 



(0,0,0), 
(1,0,0) 

(0,0,0), 
(1,1,0) 
(1,0,0), 
(0,1,0) 

(1,0,0), 
(1,0,0) 

(1,0,0), 

(-1,0,0) 

(1,0,0), 

(0,0,0) 

(1,0,0), 

(0,1,1) 



1.030 0.95 



+ 3 

-3 



1.026 0.99 



1.060 0.95 

1.062 0.91 

1.000 1.02 

1.12 0.76 +f 



+ 10 
-10 

+ 10 

-11 



1.09 0.90 



-17 
"16 



+ 2 



1.037 0.94+^ 



1.048 0.91 



+ 5 



1.050 0.96 



1.057 0.92+}° 



1.000 0.97 

1.11 0.75 

1.030 0.94 

1.08 0.91 



+11 

-12 

+ 10 

- 8 

+ 5 
-5 

+ 17 
-16 



1.073 0.93+^ 
1.069 0.89 
1.000 0.95 



1.09 0.89 



+ 8 



1.14 0.71 1^ 



1.030 0.94 



1.11 0.85 1}^ 



1.08 0.87+^° 



1.002 0.94+}^ 



1.16 0.67l^ 



1.030 0.94 



+ 5 



1.13 0.78 +}I 



Table U (continued) 
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Kn = 0.14144, K,/2 = 0.14144 



[Pmii 



0.133 0.133 



0.133 0.121 



a; 



0.133 ^ 0.125 



0.133^ 0.129 
^ ^QQ' 



(0,0,0), 
(1,0,0) 

(0,0,0), 

(1,1,0) 

(1.0.0). 

(0,1,0) 

(1,0,0), 

(1,0,0) 

(1,0,0), 
(-1,0,0) 

(1,0,0), 
(0,0,0) 

(1,0,0), 
(0,1,1) 



1.048 0.91 ^3 



1.09 0.82+^ 

1.10 0.82+1 



1.000 1.02+1^ 



1.20 0.57 



1.15 0.71 



+i'i 

-13 



1.026 1.02 ^3 



1.050 0.96 1^ 



1.074 0.88 1^° 



1.003 0.98+11 
1.15 0.69 1^? 



1.048 0.90+^ 



1.100 0.83+}^ 



1.030 1.00^3 



1.08 0.97 1^ 



1.07 0.87 1^° 
1.002 0.98+1^ 
1.16 0.67 +? 



1.048 0.91+^ 
1.111 0.80 



-5 

f 16 
-13 



1.037 0.96 ^3 



1.073 0.93^7 
1.09 0.85+^ 



1.001 0.97+13 



1.173 0.63+6 



1.048 0.92+^ 



1.13 0.76+}^ 



Table 9: Estimates of the function ^qq> as obtained from the axial form factor 
Gi. All the transitions, corresponding to initial heavy k = 0.133 and 0.129 and 
initial and final momenta up to \p\, \p'\ = \/2pmin o,re shown. Statistical errors 
in uj are in the last digit or beyond. 
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Kn = 0.14144, Kz2 = 0.14144 



bmir 



0.129 ^ 0.129 



0.129 0.121 



0.129 0.125 



0.129 ^ 0.133 



0,0,0 
1,0,0 

0,0,0 
1,1,0 
1,0,0 
0,1,0 
1,0,0 
1,0,0 

1,0,0 
-1,0,0) 

1,0,0 
0,0,0 
1,0,0 
0,1,1 



1.037 0.94 

1.073 0.93 

1.076 0.88 

1.000 1.02^ 



-3 



1.15 0.68+^ 



1.11 0.81 



-16 
-14 



1.026 1.01 : 
1.050 0.97 



1.064 0.91 +^9 
1.001 0.98+11 



1.13 0.72*° 
1.037 0.92+^ 



1.09 0.89 



1.030 0.981^ 

1.060 0.96+^ 

1.069 0.89 

1.000 0.97 

1.14 0.71 

1.037 0.92 



1.10 0.86 



-17 
-14 



1.048 0.91 
1.09 0.88 
1.09 0.84 



1.001 0.96+^3 



1.17 0.63 



1.037 0.92 



1.13 0.74 



-14 
-14 



We observe that: 



Table ^ (continued) 



• Our determinations of ^qqi{ijj) have statistical errors ranging from 3 to 15%. 
As a general rule, we note that for the same heavy-quark masses, errors 
increase with the momentum of the final baryon. Furthermore, errors are 
amplified as the light quarks approach the chiral limit. 

• The estimates of ^QQ'(t^) obtained from transitions with final momentum 
|p 'I = \/2pmin are in agreement with those of similar u) obtained from chan- 
nels with lower momentum. However, since these points are more affected 
by discretisation errors, we do not include them in our determination of the 
Isgur-Wise function. 

The estimates reported in Tables ^ and | are also plotted in Figure 0, for all 
the degenerate and non-degenerate transitions, at = K12 = 0.14144. We have 
included all the points obtained from baryons either at rest or with momenta 
|p 'I = Pmin- The interpretation of the dependence of ^qq'{uj) on the velocity 



transfer is postponed to Section |0 



We conclude this section with a comment on the determination of the suppressed 
form factors and G^. As can be observed from Figure |l|B, the signal is very 
noisy and compatible with zero. Because of this feature, which is common to 
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0.129 -> 0.129 



p,p' 


f^ll = 


0.14144, Kn = 0.14226 


/t/i = 


0.14226, Kn = 0.14226 


[Pmin] 


LO 




a; 




(0 ()) 

(1,0,0) 


1.040 


0.92+^ 


1.043 


0.92+^ 


(() ()) 
(1,1,0) 


1.078 


0.95 


1.084 


0.98 lio' 


fl,0,0), 
(0,1,0) 


1.081 


1 1 o 

0.85 ^n 


1.088 


0.80^15 


(i,0,Oj, 
(1,0,0) 


1.000 


1.02 +iJ 


1.000 


0.851^^ 


(1,0,0), 
(-1,0,0) 


1.16 


0.68+^ 


1.18 


0.72 


(1,0,0), 
(0,1,1) 


1.121 


0.77^19 


1.131 


0.59 r'l 



Table 10: Estimates of the function ^qq' as obtained from the axial form factor 
Gi. All the degenerate transitions kq — kqi — 0.129 and light masses corre- 
sponding to Kii, Ki2 = 0.14144, 0.14226 and 0.14226, 0.14226 with initial and final 
momenta up to \p\, \p'\ — -\/2pmin a'^e shown. Statistical errors in uj are in the 
last digit or beyond. 
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Figure 2: ^qq^oj), as obtained from axial (left) and vector (right) current matrix 
elements. 
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most of the momentum channels, we are unable to determine these form factors, 
within the available statistics. 

4.2.2 Analysis of the vector form factors 

The analysis of the vector form factors proceeds along the same lines as those fol- 
lowed for the axial form factors, with some significant differences. It is interesting 
to measure not only the three form factors F^^, with i = 1,2,3, separately, but 
also their sum because of the normalisation condition (p^lf). Indeed, as we will see 
below, it proves to be easier to determine the sum precisely than the individual 
form factors. 

We now study the different spinor components of the two matrices V^, W for 
different values of the index /i. As was the case for the axial current, the discussion 
is based on equations (|85D and (|86|) in Appendix ^ We find that 

1. for indices = i with i = 1,2, 3: the large components of Vj are located in 
the top-right (bottom-left) submatrix in the forward (backward) part of the 
lattice, whereas the non- vanishing components of W are located along the 
diagonal. Thus, the form factors Ff" and give separate contributions to 
different spinorial components. However, as shown in Appendix 0, spinorial 
components located in the off-diagonal submatrices are always proportional 
to the amplitudes a and f3 and thus are noisier than those located in the 
diagonal submatrices. From this feature we can already anticipate that the 
determination of Ff" is going to be less precise than that of Gf. 

2. for current indices fj, = 0: the large components of both and W are 
located in the top (bottom)-diagonal submatrix in the forward (backward) 
part of the lattice. In this case one obtains the linear combination of the 
three form factors F^ + vq F2 + v'q F^ which is approximately equal to 
Fi + F2 + -F3 at the low momenta which we are using. These channels 
have a very clean and precise signal, which is the reason for the accurate 
determination of the sum of the form factors. 

In view of the above discussion, it has proved to be advantageous to consider 
equations both with ^ = i and with /i = 0, minimizing the full covariance matrix 
with respect to the three form factors, and reconstructing their sum. Examples 
of plateaux obtained from channels with /x = and with fi = i are shown in 
Figures |I|C and |Ip, respectively. As expected, the relative statistical error is 
smaller in the first case, and we find that the stability and symmetry of the 
plateau is also better. In extracting the values of the form factors we have again 
restricted the fitting range to three time slices centered around t = 12, and we 
find very reasonable values of Xdof ■ Confirming the indications of the preliminary 
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discussion above, we measure Fi,F2 and individually with relatively large 
statistical errors, ranging from 6 to 20%, whereas the errors on their sum are 
never more than 10%. 

In order to extract the form factors from the correlation functions measured 
on the lattice it is necessary to determine the renormalisation constant Zy which 
relates the lattice and physical vector currents. This constant has previously been 
determined non-perturbatively, by studying the matrix elements of the charge 
operator between meson states, and the results are reviewed in Appendix 0. 
Noting that for degenerate transitions between baryon states at rest, Zy Vq is 
the charge operator, we find here that 

l = ^v(E^^(l))- (39) 

i 

The numerical values of Zy obtained in this way are: 



(40) 



which are in agreement with the other non-perturbative estimates presented in 
Appendix 0, obtained with a different method. Below we will use 1/ I]i FiiX) to 
normalise the lattice vector current. 
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In Tables 11, 12 and 13, we present our results for the quantity 



= THW) 

for all the combinations of quark masses and for initial and final momenta up 



to \p\, \p'\ = Pmin- The results reported in Tables |TT] and are also plotted in 
Figure 0, for all the degenerate and non-degenerate transitions, at nn = K12 = 
0.14144. The study of the dependence of ^qq/{u) on the velocity transfer is again 
postponed to Section O. 



We end this subsection with a discussion of the determination of the form factors 
Fi,F2 and F3 separately. The form factors F2 and F3 start at 0{as) in per- 
turbation theory and at next-to-leading order in the l/mq expansion and they 
are thus expected to be small relative to Fi and sensitive to the values of the 
heavy-quark masses (which in the present study vary by almost a factor two: 



m 



G [1.48,2.38]). In Figure ^ we plot our estimates of 



F2,,{u) = Fi,{u) (42) 
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Kii = 0.14144, Ki2 = 0.14144 



P,P 

[Pmii 



0.121 0.121 



0.121 ^ 0.125 



0.121 ^ 0.129 



0.121 ^ 0.133 



0,0,0), 
1,0,0) 

1,0,0), 
0,1,0) 

1,0,0), 
1,0,0) 

1,0,0), 
-1,0,0) 
1,0,0), 
0,0,0) 



1.026 1.00 



+ 3 



1.052 0.99 



1.000 1.00 



1.030 0.94 



1.057 0.87 tlo 



1.000 1.08 



1.037 0.91 



1.064 0.83 tl 



1.001 0.97 til 



1.103 0.73 



+ 9 
-9 



1.113 0.74 



1.026 0.99 



+ 8 
-7 



+ 4 

-5 



1.127 0.67 +^ 
1.026 0.98 



-7 

+ 4 
-4 



1.048 0.84 tt 

1.074 0.79 tl 

1.003 0.86 t\l 

1.145 0.61 1^ 

1.026 0.95 tl 



{Pmh 



0.125 0.125 



a; 



0.125 0.121 



a; 



0.125 ^ 0.129 



UJ 



0.125 ^ 0.133 



0,0,0), 
1,0,0) 

1,0,0), 
0,1,0) 

1,0,0), 
1,0,0) 

1,0,0), 
-1,0,0) 

1,0,0), 
0,0,0) 



1.030 0.99 +3 



1.026 0.94 tl 



1.037 0.92 tl 



1.062 0.95 



+ 9 
-9 



1.057 0.86 



+ 10 
-10 



1.000 1.00 



1.124 0.68 +1 



1.000 1.14 +!5 



14 



1.113 0.72 +^ 



1.030 0.94 tl 



1.069 0.83 
1.000 1.03 
1.138 0.66 +1 



1.030 0.97 tl 



1.048 0.85 +4 



1.080 0.78 +o 



1.002 0.91 



-11 

-13 



1.157 0.59 



1.030 0.94 +4 



Table 11: Estimates of the function ^qqi as obtained from the vector form factor 
combination J2iFi. All the transitions, corresponding to initial heavy k — 0.121 
and 0.125 and initial and final momenta up to = V^Pmm o-fe shown. 

Statistical errors in uj are in the last digit or beyond. 
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Kii = 0.14144, Ki2 = 0.14144 



P,P 

[Pmii 



0.129^ 0.129 



0.129 



y 0.121 
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^ 0.125 
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0,0,0), 
1,0,0) 
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0,1,0) 

1,0,0), 
1,0,0) 

1,0,0), 
-1,0,0) 

1,0,0), 
0,0,0) 



1.048 0.92 



1.098 0.82 



1.000 0.97 



+ 8 



-12 
-12 



1.195 0.53 



1.026 
1.074 
1.003 
1.145 
1.048 



1.00 +2 
0.86 t'l 
1 21 

-^•^-^ -15 

0.66 1^ 
0.88 +7 



1.030 
1.080 
1.002 
1.157 
1.048 



0.96 +^ 
0.84 
1.20 
0.63 
0.88 



-2 

+ 10 
- 9 

+ 13 
-15 

+ 7 

-6 

+ 6 
-6 



1.037 0.88 +3 

1.087 0.79 tl 

1.001 1.13 til 

1.173 0.58 tt 

1.048 0.87 ^5 



Table 12: Estimates of the function ^qq' as obtained from the vector form factor 
combination J2iFi. All the transitions, corresponding to initial heavy k — 0.133 
and 0.129 and initial and final momenta up to = V^Pmm o-fe shown. 

Statistical errors in uj are in the last digit or beyond. 
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0.129 0.129 



p,p' 


Kii = 0.14144, Ki2 = 0.14226 


Kii = 0.14226, Ki2 = 0.14226 


(0,0,0), 
(1,0,0) 

(1,0,0), 

(0,1,0) 

(1,0,0), 

(1,0,0) 

(1,0,0), 

(-1,0,0) 


1.040 0.95 tl 
1 nsi n R7 '^^^ 

±.UO± U.Ol -|^-|^ 

1.000 0.86 til 
1.16 0.62 +^ 


1.043 0.93 +^ 
1 nss n 77 

1.000 0.63 til 
1.18 0.64 



Table 13: Estimates of the function ^qqi as obtained from the form factor J2i Fi- 
All the degenerate transitions kq = kq' = 0.129 and light masses corresponding 
to Kii, K,i2 = 0.14144, 0.14226 and 0.14226, 0.14226 with initial and final momenta 
up to \p\, \p'\ = V^Pmin o-^c shown. Statistical errors in u are in the last digit or 
beyond. 

We observe that these form factors are small and negative as one would expect 
from Eq. (p!2D. Our measurements of 

i^iH = i^iV)|^ (43) 

for degenerate transitions are also presented in Figure ^. In order to decrease 
the number of parameters in the fit, we restricted the analysis to degenerate 
transitions, for which one can use the symmetry relation ^2(0;) = F-i{ijj). 

As anticipated above, the errors in the individual form factors are larger than 
that in their sum. To make this clear, we plot in Figure ^ the form factors Fi, F2 
and the sum, Fi + 2F2, for the degenerate channel with heavy quark hopping 
parameter kq = 0.121. It can be seen that the fluctuations in Fi and F2 partially 
compensate each other in the determination of the function C,qq' 

A more detailed discussion of the behavior of both Fi and ^2,3 and of their relation 
with the Isgur-Wise function, will be presented in Section ^.1| . 

5 The Isgur-Wise function 

In this section we study the dependence of the function Cqq'{^) on the masses of 
the initial and final heavy quarks so as to extract the corresponding Isgur-Wise 
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Figure 3: Form factors Fi and F2 = F^, as obtained from degenerate quark transi- 
tions. Different symbols correspond to different heavy-quark hopping parameters. 
The light-quark hopping parameters are always kh = K12 = 0.14144. 
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Figure 4: Example of the relative size of the normalized vector form factors, 
i^ii — 1^12 — 0.14144 and kq — kqi — 0.121. 
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function. We also study its dependence on the velocity transfer d;, and the mass 
of the light quarks. We further attempt to estimate the size of l/mq corrections. 




5.1 ^QQ'('^) as a function of to and of heavy-quark mass 

As discussed in Section H, both the axial form factor Gi and the sum of the 
vector form factors J2i Fii are protected from l/mq corrections at zero recoil (see 
Eq. ([T^)). Away from zero recoil, these quantities are no longer protected by 
Luke's theorem and they suffer l/mq corrections. We recall their expansion in 
powers of the inverse heavy-quark mass, given in Eqs. (|lT]) and (P) 

+0(l/mj(,j. (44) 

In this section we study the dependence of the form factors on uj in order to 
extract some phenomenologically interesting quantities. In particular, the slope 
of ^qq/{uj) at zero recoil can be related to the slope of the physical form factors 
(see also Section |]) through the correction coefficients given in Eqs. ( pH]) and (|12|). 
The form factors, in turn, are needed in the calculation of the decay rates and 
asymmetry parameters. 

In order to obtain reliable estimates of phenomenological quantities, we must 
learn how to extrapolate our data, obtained for initial and final heavy quarks 
with masses around that of the charm quark, to the physical b c decays. 
HQET provides us with the guide for this extrapolation, and it is important to 
understand the role of the l/mq corrections, present in the function (|4^) , and to 
check that higher-order corrections are small. 

With the aim of reducing the statistical error we exploit the relation (PD, and fit 
the vector and axial data together. This is correct up to terms of 0{l/mq) and 
two-loop perturbative corrections, which we neglect throughout this study. Near 
= 1 we expand C,qq'{uj) as a linear function of u 

^2, 



uj) = l-p'iuj-l) + 0[iiu-iy) , (45) 

and we study whether there is any dependence of the slope, p^, on the masses of 
the heavy quarks. Our results are obtained for the set of initial- and final-state 
heavy-quark masses given in Table |^, but with fixed light-quark masses around 
that of the strange {kh = k/2 = 0.14144). It can be seen from Figure ^ that there 
is no statistical evidence of a dependence of C,qq'{uj) on the heavy-quark mass. In 
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Figure 5: ^QQ'(t^) from degenerate transitions, using vector and axial current data. 
Different graphical symbols denote different heavy quark masses, which appear to 
follow the same curve. 



order to quantify this statement, we have fitted separately to the function (^51), 
each of the four data sets corresponding to degenerate transitions obtaining 
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= 0.125 
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+ 4 
-4 


at 
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confirming that no dependence on the heavy-quark mass can be detected. 

Our conclusion from this analysis is that our data for ^qq'(co') do not show any ev- 
idence of sizable l/mg corrections, within our precision. So even though ^qq'{uj) 
is not, in principle, a universal function, it appears to be a good approximation 
to the Isgur-Wise function. Thus, from now on, we will consider ^qqi{uj) to be 
our estimate of the Isgur-Wise function, ^^^"^{u). 

For the slope of the Isgur-Wise function corresponding to light-quark masses 
around that of the strange we take the result in Eq. (^61) , which will be used later 
to study the l/mq corrections to the individual vector form factors. This choice 
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is motivated by the fact that the analysis of the l/mq corrections is based on the 
data obtained from degenerate transitions. 



5.2 ^^^^(lo) as a function of the light-quark masses 

The Isgur-Wise function depends on the quantum numbers of the light degrees of 



freedom, i.e. on the so-called "brown muck" . Previous studies on the lattice El 



demonstrated that such a dependence is not negligible in the case of mesons, 
where the brown muck contains only one light quark. Thus we might expect to 
measure an even stronger dependence of the baryonic Isgur-Wise function on the 
masses of the light quarks, and we investigate whether this is the case in the 
present subsection. 

We study the dependence of ^'''^"(cj) on the light quarks by keeping the masses of 
the heavy quarks fixed at kq = k,q> = 0.129 and letting the light quarks take the 
three values listed in Table 0. By simultaneously fitting Gi{uj) and J^iFi^u) to 
Eq. {m, we find: 





= 2.4 
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= 0.14144, 


1^12 


= 0.14144 


p' 


= 2.0 


+ 5 
-5 


at 




= 0.14144, 


1^12 


= 0.14226 


p' 


= 1.7 


+ 6 
-8 


at 




= 0.14226, 


1^12 


= 0.14226 



(47) 

The comparison between these fits and the data is shown in Figures ^ for all 
three combinations of light-quark masses and also for the results extrapolated to 
the chiral limit. 

In order to obtain the slope in the chiral limit, we extrapolate the three es- 
timates of both u and ^^*^'^(a;) obtained with kq = kqi = 0.129 and with the 
combinations of light-quark masses in Eq. {^7\). We assume that u and ^'^'^^{u) 
depend linearly on the sum of the two light-quark masses, that is 

u{kq,hii,ki2) =uj{kq) + C (-^ + -^ —) (48) 

and similarly for S,^'^'^. This assumption is supported by our results for the spec- 
trum presented in Ref. 0]. The results of the extrapolation to the chiral limit, 
which are relevant for the semileptonic decay A;, AcW, are presented in Table |n 
and in Figure Our best estimate of the slope of the renormahzed Isgur-Wise 
function for the A^ baryon is 

p2 = 1.2 + 0.8-1.1 . (49) 



Using the functional form (^Sf ), we also obtain the values of u and ^'''^"(cu) for the 
semileptonic decay H^/z^, interpolating to one strange and extrapolating to 
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Figure 6: Plot of the Isgur-Wise function at fixed heavy-quark masses 
(kq = kqi = 0.129 ) and various light quark masses, down to the chiral limit. 
Both vector (crosses) and axial (boxes) determinations were used for the fits. 
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1.0 


0.85 


-25 


0.63 


+ ^ ' 

-27 




1.0429 tl 


0.93 


+ 5 
-6 


0.95 


+ 5 
-5 




1.088 +} 


0.82 


+ 18 

-15 


0.80 


+ 13 

-14 




1.175 tl 


0.75 


+ 13 
-12 


0.68 


+ 13 
-11 


Chiral/Chiral 


1.0 


1.02 


+23 
-24 


0.79 


+26 
-32 




1.0475 +^ 


0.92 


+ 7 
-7 


0.96 


+ 8 
-9 




1.097 tl 


0.90 


+22 
-20 


0.95 


+ 19 
-19 




1.194 +^ 


0.78 


+ 17 
-17 


0.74 


+ 16 
-15 


Chiral/ Strange 


1.0 


1.03 


+ 18 
-19 


0.87 


+22 
-22 




1.0439 tl 


0.93 


+ 5 

-6 


0.96 


+ 6 
-6 




1.090 +2 


0.90 


+ 17 
-16 


0.92 


+ 16 
-15 




1.179 ^4 


0.75 


+ 13 
-13 


0.70 


+ 12 
-12 



Table 14: Estimates of ^^'^^{uj) from both axial (A) and vector (V) form factors, 
at kq = kqi = 0.129 and for three combinations of the light-quark masses, as 
well as at the physical limits. The results at uj = 1 correspond to those from the 
momentum channel (1,0,0) (1,0,0). 
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one chiral quark. The numerical results are also presented in Table |T^. Our best 
estimate for the slope for the S;, baryon is 



= 1.5 tl . (50) 

In this exploratory study we have only used a very limited set of light-quark 
masses, and hence our conclusions on the dependence of the Isgur-Wise function 
on these masses are rather weak, and the extrapolation to the chiral limit is not 
very precise. This should be remedied in future simulations. 



5.3 0{l/mQ) corrections. 



In this section we attempt to extract a value of A for the A baryon, from the 
study of the l/mq corrections in the vector form factors Fi{uj) and ^2^3(0;). We 
start by recalling the relevant expressions given in Section ^ 

= N,{u)iQQ,{u;) + 0{l/ml,,,) (51) 



Fiiuj) 



the last equality being valid in the limit of equal heavy-quark masses, to which 
the present discussion is restricted. Also in this limit, we have 



1 + 



A 



1 + 



u + 1 mQ_ 

LO-l A 



u + 1 rriQ 



A 



1 



(52) 



Combining Eqs. (|5lD and (|52D and using the functions ^qqi which were determined 
from the fits to the dominant form factors Gi and Fi (see Section one can 
view the form factors Fi and F2 as functions of A alone. In fact, Eqs. (^T]) and 
(|52D, together with the coefficient functions ^(a)), evaluated at one loop order in 
perturbation theory, can be considered as our definition(s) of the binding energy. 
As mentioned in Section ^ these form factors are sufficiently sensitive to A for 
us to attempt an estimate A from the lattice data. 

To the order we are working in, the power corrections can be expressed as powers 
of the inverse quark mass or of the inverse mass of any hadron containing the 

heavy quark, provided that the same prescription is used in the evaluation of the 

^ ('5) 

coefficient functions Ci . We have decided to use the inverse quark mass, defined 
as 



m, 



Q 



M, 



A . 



(53) 



The analysis is performed using our data for degenerate initial and final heavy 
quarks (kq = kq/ = 0.121,0.125,0.129 and 0.133) for fixed light-quark masses. 
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^ iQQ' 


121 121 
0.125 0.125 
0.129 0.129 
0.133 0.133 


1 99 
1.0 0.95 
1.0 0.89 t^l 
1.0 0.81 


1 026 1 02 ^ n 
1.030 0.98 tl 
1.037 0.91 1^ 
1.048 0.83 tl 


0.121 0.121 
0.125 ^ 0.125 
0.129 0.129 
0.133 ^ 0.133 


1.050 0.95 +6 
1.060 0.91 +5 
1.073 0.83 tl 
1.093 0.71 t\ 


1.103 0.94 
1.124 0.90 
1.152 0.85 
1.195 0.74 +Jf 



Table 15: ^qqi as obtained from the form factors Fi for degenerate heavy quark 
transitions using Eqs. (|5^ and (j5^ and the value of A given in Eq. (|5^. 



f^ii = 1^12 = 0.14144. The fit to Fi{uj) is good, and we obtain 

A = 0.75 +1 GeV with yi, = 1.0. 



(54) 



The less certain fit to F2{uj) confirms this value with unexpected precision, given 
the statistical errors affecting individual points. 



A = 0.74 



tl GeV 



with 



^dof 



1.6. 



(55) 



In both Eqs. (|54D and ( ^Sf ) the first error is statistical and the second is due to the 
uncertainty in the value of the lattice spacing, see Eq. (|T^). Using these values 
of A we evaluate the coefficients A'^i^2(i^) and estimate the function ^qq>{uj) from 
the form factors Fi and F2 using Eqs. (0). These estimates are presented in 



Tables and |T^, and are compared with the functions ^qq'{uj) obtained using 
the dominant form factors, Gi and J^Fi, in Figure ^ 

Let us comment on the results for the A-baryon binding energy. 

• As expected, the form factors Fi{uj) and F2{uj) have a significant dependence 
on the masses of the heavy quarks, a dependence which can be partially 
understood in terms of the 1 / mq corrections in Eqs. (^Tj) and (|52D with the 
values of A in Eqs. ( ^if ) and (|55|) . 



We have neglected all higher-order l/mq corrections. These could alter the 
results of Eqs. (0) and (|55D by as much as 30%. 



The separation of l/rriQ corrections from discretisation errors is difficult. 
Although we have used an improved action to reduce these errors, and have 
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0.121 0.121 
0.125 0.125 
0.129 0.129 
0.133 ^ 0.133 


1.0 0.96 
1.0 0.79 ttl 
1.0 0.61 t% 
1.0 0.42 


1.026 1.10 +17 
1.030 0.92 +1^ 
1.037 0.75 +^5 
1.048 0.58 


0.121 0.121 
0.125 0.125 
0.129 0.129 
0.133 ^ 0.133 


1.050 1.15 
1.060 0.91 +3^ 
1.073 0.68 +28 
1.094 0.47 +22 


1.103 1.87 
1.124 1.79 til 
1.152 1.63 
1.195 1.38 til 



Table 16: ^qq' as obtained from the form factors F2 = F3 for degenerate heavy 
quark transitions using Eqs. and ([5^ and the value of A given in Eg. (|5^. 



fqq, from 




1.0 1.1 1.2 1.3 1.0 1.1 1.2 1.3 



Figure 7: ^qq'(w) as obtained from the form factors Fi and F2 = F^,. The points 
are compared with the best fit of C,qqi{uj) obtained previously from the dominant 
form factors, Gi and J2 Fi. 
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normalized the form factors by Zy defined by Eq. (p9D, the stabihty of the 
results as the lattice spacing and action are varied should be checked in 
future simulations. The consistency of the two results (|5^) and (|55D may 
be evidence that discretisation errors are relatively small. 

Although we have only performed this study with light quarks with masses 
corresponding to = K12 = 0.14144, we can nevertheless calculate Aa and 
Ah by using 

A-Aa = MBaryon - Ma + C»( — ) 

and (56) 

A - Ah = MBaryon " Mh + 0{ — ) , 

where A and M-q^^.^^^ are the binding energy and mass of the baryons 
which we study in our computation. In this way we find: 

Aa = 0.37 til GeV and Ah = 0.50 GeV , (57) 

where statistical and systematic errors have been combined in quadrature. 
The value of Aa would correspond to a mesonic binding energy of Ames — 
80 ± 120 MeV, in agreement, within the large error, with the result quoted 



in Ref. 0. The central values of the quark masses obtained using Eq. (|57|) 
are different from those used throughout this paper which were obtained 
from Eq. (p^B]). We have checked that this difference does not noticeably 
affect the correction coefficients CI , and hence our estimates of physical 
quantities. 

In addition to the systematic uncertainties arising from lattice artefacts, it 
must be remembered that the computation of power corrections in general 
is a complicated subject (for a recent review see |]22|). In the present case, 
following standard procedure, we are trying to quantify l/mq corrections 
when we are ignorant of the 0{a'j.) terms in the perturbation series for 
the coefficient functions. Although, within our errors, we have found no 
inconsistencies, we cannot be sure that the values of A would not change 
significantly if higher-order perturbative terms were included, or that the 
values would agree with other definitions of A. 



6 Phenomenological implications 

In this section we use the results for the baryonic Isgur-Wise function and for the 
baryon binding energy A computed above, to obtain the physical form factors 
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corresponding to the semileptonic decays: 

Afc ^ Ac + ev and S;, ^ Sc + ez/ . (58) 

We also derive the expressions for the decay rates near zero recoil and use them 
to make quantitative predictions. We will discuss in detail the decay of the A;,; 
that of the is very similar. 



6.1 Physical form factors 

The physical form factors can be reconstructed from the computed Isgur-Wise 
function ^''''"(o;), using the relations in Eq. (|ll]). In the following we will use 
the result, discussed in Section ^7\[ that the function ^qq'(co') is effectively in- 
dependent of the mass of the heavy quark. Therefore, the form factors de- 
pend on the masses of the heavy quarks only through the correction coefficients 
Nf\uj) = N-^\uj,mQ,mQ'), which are calculated using expressions (JT^). The 

short-distance coefficients, cl^\u!,mb,mc), were computed for Aqcd = 250 MeV, 
Uf = 4 and by fixing the quark masses to the values 

rub = 4.8 GeV, = 1.45 GeV. (59) 

The factors CI depend on the quark masses very mildly. On the other hand, the 
correction coefficients N- are very sensitive to the quark masses, because they 
contain l/mg corrections. In this case, we have expressed nib and rric as 

mb = Ma, - A and = Ma, - A (60) 

with Ma, = 5.64 GeV and Ma, = 2.285 GeV, respectively. We report our 

(5) 

estimates of the correction coefficients N- for various values of oo, for the A and 
S decays, in Table 

It is convenient to expand the physical form factors in u; — 1, near zero recoil: 



Fi{u;,mb,mc) =r]Y - p]^{uj -1), Gi{uj,mb,mc) = r]f ~ pf{uj - 1); (61) 
ler 

Nl and to the slope of the Isgur-Wise function by 



where the normalisations riY'^ and the new slopes are related to the coefficients 



r]Y = N,{l,mb,m,), r]f = Nf{l,mb,m,), (62) 

dNi{u, nib, 



Pi = p Ni{l,nib,nic) - 



(63) 

UJ=l 



and 

(64) 



~A 2 Ar5/i N dNf{u;,nib,nic 
Pi = P N^{l,mb,ni^) — 



1.^=1 
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decay 




A^i N2 N3 Nf Nl Nl 


Ac 


1.0 
1.1 
1.2 
1.3 


128 """^ -0 19 +1 -0 06 99 -0 24 +5 09 
1.25 t\ -0.18 +4 -0.06 +1 0.97 -0.23 t\ 0.08 +2 
1.21 t\ -0.17 1;] -0.05 t\ 0.95 -0.21 +^ 0.08 +^ 
1.15 1^ -0.16 +3 -0.05 t\ 0.91 -0.19 0.07 




1.0 
1.1 
1.2 
1.3 


1.32 +^ -0.23 t\ -0.07 +} 0.99 -0.28 t\ 0.10 +2 
1.30 +^ -0.22 t\ -0.07 +} 0.97 -0.26 t\ 0.10 +2 
1.26 t\ -0.20 -0.06 1} 0.95 -0.24 t\ 0.09 +2 
1.19 +5 -0.18 +3 -0.06 t\ 0.91 -0.22 t\ 0.08 t\ 



Table 17: Correction factors needed to relate the form factors at the physical limit 
with the Isgur-Wise function. 



decay 








F2 




F3 






G2 




G3 




Ac Ac 


p 


1.8 1 


0.9 
1.5 


-0.4 


+ 2 

-1 


-0.10 


+ 7 

-4 


1 Q +0.8 
i.O -1.2 


-0.4 


+ 3 
-2 


0.16 :t 


- 6 
10 


1 — :c 1 — 'c 


p 


2.3 1 


0.9 
1.4 


-0.5 


+ 2 

-2 


-0.15 


+ 7 

-4 


1 6 

-1.0 


-0.6 


+ 3 
-2 


0.22 1 


- 7 
10 



Table 18: Slope of the physical form factors near zero recoil. 



Our results for pj'^ are presented in Table 0. We observe that this procedure 
has the effect of taking us back to the form factors for physical 6 — c decays from 
the Isgur-Wise function, which, in turn, was determined by dividing the lattice 

data for the form factors (for unphysical quark masses) by the coefficients N^^\ 

(5) ^ 
Clearly, most of the uncertainty in the factors N- , due to their dependence on A 

and the quark masses, is now cancelled, as it should be for any physical quantity. 

To quantify this statement, we report the result of the following exercise. We 

have measured the ratio of the slopes of the form factor Gi{uj), at the chiral 

limit, letting A vary from 200 to 600 MeV, and changing the coefficients A^j 

accordingly. We obtain 



(A = 200) 
pf(A = 600) 



1.00 I'l , (65) 



where the relatively large error is largely due to the extrapolation to the chiral 
limit. 
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6.2 Decay rates 



Following Refs. |23] and |24], we define the helicity amplitudes, in terms of the 



physical form factors (0), in the velocity basis: 



,y,A ,/2Ma,Ma^(^t1; 

X ((Ma, ± MaJFi^'^ ± Ma,(cu ± l)^^'''^ ± Ma, (a; ± 1)^3'''^) 



i^r-, = , ' ° X (66) 



HY'f = -2JMa,Ma,(cut1)Fi^"^ 
where, for brevity, we set 

Fr = F,{u;), Ff = G,{u). (67) 

and where the upper sign corresponds to V and the lower one to A. 

The helicity amplitudes H^'J]^^ carry information about the helicity of the current 
(AvK = for a longitudinally polarised W and Xw = ±1 for transversely polarised 
one), and of the daughter baryon Ac (Ac = ±1/2). The missing amplitudes can 
be computed by means of the relations 

hY a =±hYx ■ (68) 
For convenience, we also define: 

Hx,,xw = Hl^,\w + ^t^w ■ (69) 

Differential decay rates can then be evaluated: 







duj 


(27r)3 






duj 


(27r)3 



1'^'' UMl 



where and Tl are the contributions to the rate from transversely and longi- 
tudinally polarised W^s respectively, and whose sum is 



du {271^^ ''^ 12Ma, ^ ^ 

X [\Hi/2,lf + |if_l/2,-l|^ + \Hi/2,of + \H_i/2fi\^ 

As can be seen from Eqs. (|70|) and ([7ll), these quantities can be estimated, near 
zero recoil, using our results for the form factors, in a model-independent way. 
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In complete analogy with what is done for B D^*^iu decays, we define a form 
factor, B{uj), which reduces, in the heavy-quark limit, to the Isgur-Wise function, 
^'''^'^(ci;), defined in Eq. (|^). In terms of this form factor, the rate of Eq. (|7TD is: 



„2 



uj+l\ 1 + r2 -2r(2cj+ l)/3 



2 J 1-r 



|2 



/ /cj-l\ l+r2-2r(2LU-l)/3\ ,,2 

with r = M\^/M\^. Near u; = 1, the form factor can be expanded as 

B{lo) = G,(l){l-pl{u-l) + {{to - If) } (73) 

and the results of Tables [l^ and [l^ can be used to determine G'i(l) and p|. We 
find, combining errors in quadrature, a slope 

pI = 1.1 ± 1.0 (74) 

for Ab — > Acii' decays and 

p| = 1.4±0.8 (75) 

for "Eciv decays. These values are the ones that should be compared to 

the slopes obtained by performing fits to experimental results for dV/duo versus 
a;, for a; near 1, when such results become available. The results of Eqs. ( |7^ ) 
and ([75|) can also be compared to the slopes we found for the corresponding 
Isgur-Wise functions (Eqs. (|^) and (|50D). These two sets of slopes are virtually 
indistinguishable, especially given the size of our present errors. 

For both decays we find Gi{l) = 0.99 which is just 1 + 5f (l)^^, as it should be at 
the level of precision at which we are working (see Eqs. (H) and (p!^)). In princi- 
ple, though, Gi(l) receives also I/ttlq corrections and higher-order perturbative 
corrections, both of which are beyond the precision reached in the present paper. 



We now turn to integrated rates. The physical limit for lu extends up to cij 1.43, 
which is beyond the range of velocity transfer accessible to us [u G [1.0, 1.2]) in 
the present simulation. We thus define the partially-integrated decay rate, 

rr'^*(t^max) = ^ rfcj^ (76) 

as a function of the upper limit of integration, for each of the rates i = T,L and 
T + L. 

In Table we present our results for the quantities rf'^'^*(ci;inax) for several values 
of i-i^max- For the case of the A^, the A;, and the we have used the experimental 
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decay 
A. 



ppart 



1.1 

0.57 



+ 9 

-7 



1.15 
0.98 + 



25 
18 



1.20 
1.4 



1.25 



+ 5 
-4 



1. 



+ 9 

- 7 



1.30 
2.2 



+ 1.4 

-1.2 



^max 
ppart 



1.1 

0.66 



+ 7 

-8 



1.15 
1.1 



+ 2 
-2 



1.20 

1.6 + 



1.25 
1.9 +1 



1.30 
2.2 



+ 1.2 
-1.3 



decay 
A,, ^ A, 



^max 
ppart 



1.1 

0.23 " 



1.15 
0.44 + 



1.20 
0.71 



+ 17 
-13 



1.25 
1.0 tl 



1.30 
1.4 + 



^max 
ppart 
J- L 



1.1 

0.28 



1.15 
tl 0.54 tl 



1.20 



0.86 



-14 
-16 



1.25 
1.2 +: 



1.30 
1.7 + 



decay 

Af, Ar 



^max 
ppart 



1.1 

0.34 " 



1.15 
0.53 + 



1.20 
0.7 + 



1.25 
0.8 



1.30 



^max 
ppart 



1.1 
0.38 



+ 5 

-5 



1.15 

0.58 + 



1.20 
0.7 



1.25 



+ 3 

-3 



1.30 



Table 19: Partial decay rates, in units of iKhp 10"*^^ s^^, for the A and H semilep- 
tonic decays for various values of uj^^^. The transverse decay rate is very sensitive 
to quadratic terms in {uj — 1), and the predictions foruj > 1.2 are no longer reli- 
able. 

values for the masses, whereas, for the Eh, which is as yet undiscovered, we have 
used the value computed in our previous paper on heavy baryon spectroscopy : 

Ms, = 2.47 GeV [Exp] Mh, = 5.76 GeV [Latt]. (77) 

At present, a direct comparison of our results with experiments is not possible. 
Although the semileptonic decay of A^ has been observed by various experiments 
[0], a measurement of the decay rate is not yet available. The problem of deter- 
mining the rate of the Ab and Eb semileptonic decays has been addressed theo- 
retically, making use of different models and approaches (e.g. infinite momentum 
frame (IMF), dipole form factors in Ref. and quark model in Ref. p4|), by 



several authors. Their predictions for the rate, integrated up to the end-point, 
are reported in Figure ||, and compared with the function r(ct;max)- To evaluate 
this function we have assumed \Vcb\ = 0.044. 

Finally, we note that many other interesting quantities could be computed and 
confronted with future experiments, e.g. asymmetry parameters (see for example 
[25]) and the ratio of the longitudinal to transverse rates. However, all these 
quantities are sensitive to either I/tuq effects or to terms proportional to (cu — 1)^. 
Both of these effects are beyond the precision reached in the present study. 
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Figure 8: Partially-integrated decay rate for the process — >• + W, in units of 
10^° s~^, as a function of the limit of integration uj^ax- A comparison with several 
model estimates is shown at the end-point ^|/. The solid curve corresponds 
to our central values for p and the dotted curves to the errors on these values. 
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Conclusions 



During the past few years lattice simulations have been applied very success- 
fully to weak decays of heavy mesons, and in this work we have extended these 
techniques to heavy baryons. We have presented an extensive lattice study of 
the semileptonic A;, — > Ac + IP and Hf, — > Sc + IP decays, resulting in predic- 
tions for the decay rates, for values of the velocity transfer up to about uj = 1.2. 
We have developed the formalism necessary for extracting the decay amplitudes, 
demonstrated the feasibility of obtaining phenomenologically interesting results, 
and presented the first set of predictions. We anticipate that the apphcation of 
lattice QCD to studies of the forthcoming experimental data on heavy-baryon 
decays will be an active area of phenomenology during the coming years. 

HQET is an important tool in the application of lattice QCD to weak decays of 
heavy baryons, as it is also for mesons. In order to keep lattice artefacts reason- 
ably small, we are forced to perform the computations with heavy quarks with 
masses not much larger than that of the charm quark, and then to extrapolate 
the results to physical b ^ c decays. HQET provides us with a sound theoretical 
formalism for performing this extrapolation. 

Perhaps the weakest feature of our study was the inability to determine the 
behaviour of the decay amplitudes with the mass of the light quarks with sufficient 
precision. We did perform the computations for three combinations of masses for 
the two light quarks, which allowed us to attempt an extrapolation to the physical 
A (with two almost massless valence quarks) and S (with one strange and one 
almost massless quark) baryons. However, in performing these extrapolations 
the statistical errors are amplified, and one of the priorities of a future simulation 
should be to generate substantial datasets for a larger set of light-quark masses. 
For example, it will be very interesting to check whether the slope of the Isgur- 
Wise function decreases as the masses of the fight quarks are reduced, as expected 
theoretically. In the present simulation we have a hint that this is the case, but 
the statistical errors are too large to draw a definitive conclusion. 
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A The use of extended interpolating operators 

In order to enhance the signal for the baryonic correlation functions, the light 
and heavy quark propagators have been computed using the Jacobi smearing 
method |jl2|. Since smearing is not a Lorentz-invariant operation, it alters the 
transformation properties of the correlation functions, particularly at non-zero 
momentum. In this Appendix we present the formalism required to extract the 
form factors from two- and three-point correlation functions computed using ex- 
tended (smeared) interpolating operators for the baryons. 

Consider the local operator 0{x) defined in Eq. (|2^) , 

0,ix) = ea,c {lfix)C^,l',ix))Q;ix) , (78) 

where p is a spinor index. Here we have suppressed the index Q in labelling the 
operator. O has non-zero overlap with spin 1/2 states, such the A— baryon: 

{0\O,mp,r) = Zu'-;\p). (79) 

where r is the polarisation index.! The ket in Eq. ([7^) represents a heavy A state 
(e.g. Af, or Ac). The amplitude Z is a Lorentz scalar. 

The smeared baryonic operator can be written as 



0;{x,t) = e,,, 5: /(|y-f|)/(|z-x|)/(|u;-f|) 



p 

y,z,w 



X 



X {lf{y,t)C^,ll{z,t))Q;{td,t). (80) 

Because the smearing is performed only in the spatial directions, Lorentz sym- 
metry is lost and only spatial translations, rotations, parity and time reversal 
survive. Therefore, the overlap of the operator Op{x, t) with the state \p, r) is 
given by the more general expression 

{o\o;mp,r) = [{z,m + z2{\p\hoy^\p) 
(o|o;(o)|p,r) = [v^^\p}{z,m - z^mio) 

where the amplitudes Zi and Z2 may depend on the magnitude of the three- 
momentum of the state \p,r), in accord with the restricted symmetries of the 
system. 



^Our spinors are normalized such that — = [E/m)5^ 
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A.l Smeared two-point functions 



We now study the consequences of the above discussion in the case of smeared 
source and sink (SS) two-point functions. Using Eq. (pTD, we can derive the 
general expression for the two-point function at large values of t and (T — t) i 



Gf^{t,p) = E e~'^^-^0,{x, t)a(0, 0)) 



f82) 



-ip-x 



e-^(^'^*+^'^-"> I Op(0 , 0) I g, r ) (g, r I a (0, 0) 1 0) ■ 



per 



x«M(p-)(Zi(b1) + Z2(b1)7o; 
xt;M(-p-)(Zi(b1)-Z2(b1)7o' 



We find it convenient to write the spin matrix G^^ in terms of the parameters 



Zs — Zi + Z2, 
rather than Zi and Z2: 



a = {Zi- Z2)/ {Zi + Z2 
E + m — a"^ {E — m 



— e 



^E(p){T~t) 



AE 

E + m + a^{E — m) 
AE 

E + m — a'^{E — m 
AE 

E + m + a'^{E — m) 



1+ 



2a ^ ^ 



(83) 



AE 



2a ^ ^ 



Local operators, and full 4-dimensional cubic symmetry, correspond to the case 
a = 1, i.e. Z2 = 0. 



A. 2 Smeared three-point functions for the study of the 
semileptonic decay of the A/,. 

We now present the expressions for the smeared-smeared three-point functions 
from which the form factors are extracted. The V — A weak current is of course 

^Note that we are using anti-periodic boundary conditions in time. 
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a local operator; it is the interpolating operators for the baryons in Eq. (^) which 
are now smeared. In the forward half of the lattice (F), for large t, so that only 
the lightest state contributes, we obtain 



C^{t,,ty) = K{t,,ty) {[e + M + a{M - E) 



1 + 



^5) 



+ 



7o + 



(1 - a)p- 770 



X 



X 



E + M -a{M-E) 
- (1 + a)p- 7 I X 

- '{G^i^h, + G^{u)v'^ + Gl{u)v,) 75] } X 
E' + M' + /3{M' - E')] 1+\e' + M' - p{M' - E') 



{ 



7o 



(l + /3)p'-7]} 



- (1 - (3)p' ■ 77o 
whereas in the backward half (-B), for large T — t, 

= K{T-t,,T-ty){[E + M + aiM-E) 



1 



(86) 



+ 



7o - 



X 



X 



E + M-a{M -E) 
- (1 + a)p- 7 } X 



(1 - a)p- 77o 



+ 



+ 



(Gf (a;)7, + + Gl{u)v,) 75] } x 



{ [e' + M' + (3{M' - E')] 1-[e' + M' - (3{M' - E')] 70 + 



+ 



(1 - ■ 770] - [{1 + /3)p' ■ i]} 



where v = {E, —p)/M (similarly for v') and where in both these equations the 
normalisation factor is given by: 



16EE' 



Since tx = 24 = T/2 in our simulations, the dependence on is the same in both 
halves of the lattice. In Eqs. (|^) and (|86D a and /5 are the functions defined 
in (P^D, for the final and initial particle respectively. They are identical only 
for degenerate transitions when \p\ = \p'\. For local operators, a = (3 = 1 and 
expressions (^3|-^) reduce to Eq. (P?]). 

Expressions (^) and (^) are used in Section ^ to extract the values of the form 
factors from the correlation functions. Further clarification of the procedures 
which are used is pressented in the following appendix. 
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B On the analysis of three-point correlators 



In this appendix we discuss in greater detail the analysis procedure used to ex- 
tract the six form factors from the correlation functions. Such an analysis is 
complicated both by the very large number of non-vanishing components of the 
three-point correlation function, which has both spinorial and Lorentz indices, 
and by the effect of the smearing on the baryonic operators discussed in the 
previous appendix. It proved to be convenient to restrict the analysis to those 
components which are proportional to large, and hence precisely measured, kine- 
matical coefficients, as we now explain. 

For illustration we consider here the forward half of the lattice; the extension to 
the backward half is straightforward. We consider two typical cases. 



Example 1: Coefficient of the form factor for /i = 0. 



We rewrite the expression in Eq. (|85|) , which is a 4 x 4 matrix in spinor 
space, in terms of 2 x 2 matrices i: 



Vectoro 



= 4[(E + M)(E' + M')+p-p'](I + 7o) 

+ 4a/5[-(^-M)(^'-M')-p-p"](I-7o) 

+ Aip,p'^a'^{l + 7o) - Aa^ipip'^a'^it - 70) 

+ 4a[(M-E)p,-(E' + M')p.]7.(I + 7o) 

- 4/3 [(M + E)p,-(E'-M')p.] 7^(1 -7o), 



^9) 



where Vectoro denotes the time component of the three point correlation 
function with the vector current. So the matrix structure shared by the 
coefficients of the vector form factors is of the form: 



Vectoro 



OC 



V 



a 



/3 



a/3 



(90) 



where each of the submatrices is a 2 x 2 matrix. 



^The following representation for the gamma matrices is used here: 

70 = 1 n ^ l^=[ ^„ n' ) with i = 1,2, 3. 
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• Example 2: Coefficient of the form factor Gf for /i = /c = 1, 2, 3. 
A similar manipulation of Eqs. ( pSf ) for the axial current leads to 




= {4[(M + E)(M' + E')]7fe(I + 7o) (91) 



+ 4«/?[(M-E)(M'-E')]7fc(I-7o) 
+ 4a [(M - E)p[^kli - (M' + E')paak] (I + 7o) 
+ 4/3 [(M + E)p:7fc7, + {E' - M')pazlk] (I - 7o) 
+ 4 PiPj'jaklj (I + 7o)+4a/3 PiPj'jaklj (I-7o)}75, 



which also has the structure of (^). 
The relative sizes of the components of the matrices ( ]89| ) and ( ]9T| ) are as follows: 



• the submatrices proportional to the identity are proportional to MM' (plus 
small terms quadratic in the momenta), and are therefore large and precisely 
determined. 

• the submatrices proportional to a and f3 are proportional to terms like Mp'^ 
or M'pi and are therefore medium-sized. Furthermore they have an addi- 
tional statistical uncertainty due to the presence of the amplitude factors a 
or p. 

• the submatrices proportional to a/3 are proportional to [E' — M'){E — M) or 
other terms which are quadratic in the momenta, and are therefore small. 

Finally, we note that, as a general feature, the vector current will give large 
contributions for /i = and the axial current when /i = i, z = 1,2,3. 

C Determination of Zy 

In this appendix we discuss the normalisation of the lattice vector current, used for 
extracting the vector form factors. The normalisation factor, Zy has previously 
been determined from matrix elements between heavy meson states, using the 
same configurations, and for the same values of the quark masses as those used 
in the present study and we present the results below. We start however 
with a determination of Zy from matrix elements between heavy baryon states. 

Zy can be measured from the correlators (0), computed for degenerate initial 
and final heavy quarks {Q = Q'). For instance, with Dirac index i = 1,2, we 
define 



Zy{t) = /C 



C^{p,T/2) 



1 



(92) 



C{p,q = 0,T/2,t) 



At=0 



FHl) + Fi^{l) + Fi{l) 
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with 

+ M) 2E 

2E E + M-a\M-E) ^^^^ 

and where was defined in Eq. (^1]) and t is taken in (0,T/2). Zy is then 
obtained by fitting Zy{t) to a constant in the plateau region. 

We have computed the necessary two- and three-point correlators for two different 
forward channels: p = and p = (p^in) O5 0). From the correlation functions with 
the baryon at rest we obtain: 



(94) 



The statistical errors are very small (although not as small as for meson states), 
which is not unexpected since we are studying the effects of the charge operator. 

For Zy measured from correlation functions with p = (Pmiru0,0), the statistical 



Zy 


= 0.898 


+ 3 
-2 


at 




= K,12 


= 0.14144, 


f^Q 


= 0.133 


Zy 


= 0.924 


+ 2 

-2 


at 




= 1^12 


= 0.14144, 


I^Q 


= 0.129 


Zy 


= 0.948 


+ 2 

-2 


at 




= 1^12 


= 0.14144, 


I^Q 


= 0.125 


Zy 


= 0.970 


+ 1 
-2 


at 


K-ll 


= K,12 


= 0.14144, 




= 0.121 



errors are too large for us to make any comparison with the results in (p^). In 
this case we find: 



(95) 



We now compare the results for Zy obtained between baryonic states (Eq. (|9^) ), 
and those obtained in Ref. |^ from matrix elements between pseodoscalar states 



Zy 


= 0.88 


+ 15 
-10 


at 


/«zi 


= 1^12 


= 0.14144, kq 


= 0.133 


Zy 


= 0.90 


+ 14 
-11 


at 


/«/i 


= 1^12 


= 0.14144, Kq 


= 0.129 


Zy 


= 0.94 


+ 14 
-11 


at 


i^ii 


= 1^12 


= 0.14144, Kq 


= 0.125 


Zy 


= 0.97 


+ 15 
-12 


at 


i^ii = 


-- 1^12 = 


-- 0.14144, Kq = 


= 0.121 . 



using the same configurations and quark masses: 



Zy 


= 0.8913 


+ 2 

-1 


at 




= 0.14144, KQ 


= 0.133 


Zy 


= 0.9177 


+ 3 
-2 


at 


1^1 


= 0.14144, KQ 


= 0.129 


Zy 


= 0.9428 


+ 4 
-2 


at 




= 0.14144, KQ 


= 0.125 


Zy 


= 0.9659 


+ 6 

-3 


at 


f^i = 


-- 0.14144, KQ = 


= 0.121 . 



(96) 



The agreement of the results obtained for Zy using mesonic and baryonic corre- 
lation functions to within less than 1% is reassuring. 



The variation of the values of Zy with the mass of the heavy quark in ( |9^ and 
(|96|) is an effect of the discretisation errors, due to the fact that the quark Q 
is heavy. To see this more clearly we compare the results with those obtained 
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between light pseudoscalar mesons (with degenerate valence quarks): 



The results in ( pT]) were obtained on a subset of 10 gluon configurations |20|. The 
dependence on the masses of the light quarks is seen to be very mild, and the 
reuslts are in consistent with the expectations from one-loop perturbation theory 
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Zv = l- 0.10^^ + C»(/) - 0.83 ai 13 = 6.2 (98) 

when evaluated using the boosted value of the coupling constant, obtained from 
the mean field resummation of tadpole diagrams p7[ . 



The results for Zy in Eq. (^) and (^), obtained using heavy baryon and meson 
states, differ from those obtained with light mesons, (0), by about 10-20% for 
the range of quark masses used in our simulations. This difference is a good indi- 
cation of the size of mass-dependent discretisation errors in our calculation; it is 
consistent with our expectation that they should be of O (asamq) and O (^a'^mg^ . 
However, as explained in Section ^, the errors in the computed form factors are 
expected to be considerably smaller, because we normalise all of them by Zy also 
determined between heavy baryon states. 
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